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! ABSTRACT 

Aims. We examine the role of rotation on the evolution and chemical yields of very metal-poor stars. 

Methods. The models include the same physics, which was applied successfully at the solar Z and for the SMC, in particular, 
shear diffusion, meridional circulation, horizontal turbulence, and rotationally enhanced mass loss. 

Results. Models of very low Z experience a much stronger internal mixing in all phases than at solar Z. Also, rotating models 

■ at very low Z, contrary to the usual considerations, show a large mass loss, which mainly results from the efficient mixing of 
' the products of the 3a reaction into the H-burning shell. This mixing allows convective dredge-up to enrich the stellar surface 

in heavy elements during the red supergiant phase, which in turn favours a large loss of mass by stellar winds, especially as 
• rotation also increases the duration of this phase. On the whole, the low Z stars may lose about half of their mass. Massive stars 

initially rotating at half of their critical velocity are likely to avoid the pair-instability supernova. The chemical composition 
, of the rotationally enhanced winds of very low Z stars show large CNO enhancements by factors of 10'^ to 10^, together with 
' large excesses of ^^C and ^'^O and moderate amounts of Na and Al. The excesses of primary N are particularly striking. When 

these ejecta from the rotationally enhanced winds are diluted with the supernova ejecta from the corresponding CO cores, we 

find [C/Fe], [N/Fe],[0/Fe] abundance ratios that are very similar to those observed in the C-rich, extremely metal-poor stars 
' (CEMP). We show that rotating AGB stars and rotating massive stars have about the same effects on the CNO enhancements. 
. Abundances of s-process elements and the ^^C/^^C ratio could help us to distinguish between contributions from AGB and 

■ massive stars. 
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1. Introduction 

Interest in the evolution of extremely metal-poor stars has 
been stimulated recently by at least two types of observ- 
ing programmes. First, the detection of very far galaxies 
at redshifts well beyond 6 (see e.g. Pello et al. l2005|l opens 
the way to detection of galaxies whose colours will be 
dominated by extremely metal-poor stars fSchaerer i2002il 
I2003|l . Second, as a complement to the observation of 
the deep Universe, the detection of nearby, very metal- 
poor halo stars provides very interesting clues to the 
early chemical evolution of our Galaxy (Beers et al. lT^ 
Beers [TUM Christ lieb et al. IMJl Bessell et al. IMHI 
Cayrel et al. HUM Spite et al BUD5i|l . These works have 
shown the following very interesting results: 

— The measured abundances of many elements at very 
low metallicity present a very small scatter ( Cayrel et 
al. \MJ^ . At first sight this appears difficult to under- 
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stand. Indeed, in the early Universe, stars are believed 
to form from the ejecta of a small number of super- 
novae (may be only one). For instance the Argast et 
al. models f^MMVM^ predict that for [Fe/H] < -3, 
type II supernovae enrich the interstellar medium only 
locally. Since the chemical composition of supernova 
ejecta may differ a lot from case to case, large scatter 
of the abundances is expected at very low metallicity. 
For most of the elements, however, this strong scat- 
ter is not observed, at least down to a metallicity of 
[Fe/H]- -4.0 (Cayrel et al.'5DD3). This might indicate 
that, already at this low metallicity, stars are formed 
from a well-mixed reservoir composed of ejecta from 
stars of different initial masses. 
— These observations also show that there is no sign 
of enrichments by pair-instability supernovae, at least 
down to a metallicity of [Fe/H] equal to -4- Let us 
recall that these supernovae have very massive stars 
as progenitors, with initial masses between approxi- 
mately 140 and 260 Mq (Barkat et al. [TM7I Bond et 



2 



G. Meynet, S. Ekstrom, & A. Maeder: The early star generations 



al. 119841 Heger & Wooslev l!jUU2|l . Such massive stars 
are believed to form only in very metal-poor environ- 
ments. At the end of their lifetime, they are completely 
destroyed when they explode as a pair-instability su- 
pernova. In this way they may strongly contribute to 
the enrichment of the primordial interstellar medium 
in heavy elements. The composition of the ejecta of 
pair-instability supernovae is characterised by a well 
marked odd-even effect and a strong zinc deficiency. 
These two features are not observed in the abundance 
pattern of very metal-poor halo stars. Does this mean 
that at [Fe/H] equal to -4, pair-instability supernovae 
no longer dominate the chemical evolution of galaxies 
or that such stars are not formed ? If formed, could 
these stars skip the pair instability or have different 
nucleosynthetic outputs from those currently predicted 
by theoretical models ? 

— The N/0 ratios observed at the surface of halo stars by 
Israelian et al. i2004)) and Spite et al. 1^005) indicate 
that important amounts of primary nitrogen should be 
produced by very metal-poor massive stars ( Chiappini 
et al. \2()()5]) . The physical conditions for such impor- 
tant productions of primary nitrogen by very metal- 
poor massive stars remain to be found. 

— // most stars at a given [Fe/H] present a great ho- 
mogeneity in composition, a small group, comprising 
about 20 - 25% of the stars with [Fe/H] below -2.5, 
show very large enrichments in carbon. These stars 
are known as C-rich extremely metal-poor (CEMP) 
stars. The observed [C/Fe] ratios are between ~2 and 
4, showing a large scatter. Other elements, such as ni- 
trogen and oxygen (at least in the few cases where 
the abundance of this element could be measured), 
are also highly enhanced. Interestingly, the two most 
metal-poor stars known up to now, the Christlieb star 
or HE 0107-5240, a halo giant with [Fe/H]=-5.3, and 
the subgiant or main-sequence star HE 1327-2326 with 
[Fe/H]=-5.4 (Frebel et al. EUUSIl belong in this cate- 
gory. To explain such high and scattered CNO abun- 
dances, obviously a special process has to be invoked 
(see Sect. 6 below). 

The results outlined above clearly indicate that new sce- 
narios for the formation and evolution of massive stars at 
very low Z need to be explored. 

Among the physical ingredients that could open new 
evolutionary paths in very metal-poor environments, rota- 
tion certainly appears a very interesting possibility. First, 
for metallicities Z between 0.004 and 0.040, the inclusion 
of rotation improves the agreement between the models 
and observations in many respects by allowing us to re- 
produce the observed surface abundances (Heger & Langer 
I^Unni Meynet & Maeder^^, the ratio of blue-to-red su- 
pergiants at low metallicity (Maeder & Mevnet r2001|l . the 
variation with the metallicity of the WR/0 ratios and of 
the numbers of type Ibc to type II supernovae (Meynet 
& Maeder 12005(1 . Most likely, stars are also rotating 
at very low metallicity, and one can hope that the same 



physical model assumptions that improve the physical de- 
scription of stars dX Z > 0.004 would also apply to the 
very low metallicity domain. 

Second, if the effects of rotation are already quite sig- 
nificant at high metallicity, one expects that they are even 
more important at lower metallicity. For instance, it was 
shown in previous works that the chemical mixing be- 
comes more efficient for lower metallicity for a given ini- 
tial mass and velocity (Maeder & Meynet 120011 Meynet 
& Maeder ,2002). This comes from the fact that the gradi- 
ents of f2 are much steeper in the lower metallicity models, 
so they trigger more efficient shear mixing. The gradients 
are steeper because less angular momentum is transported 
outwards by the meridional currents, whose velocity scales 
as the inverse of the density in the outer layers (see the 
Gratton-Opick term in the expression for the meridional 
velocity in Maeder & Zahn lTMHji . 

Third, rotation can induce mass loss in two ways. The 
first way, paradoxically, is linked to the fact that very 
metal-poor stars are believed to lose little mass by ra- 
diatively driven stellar winds. Indeed, in the radiatively 
driven wind theory, the mass loss scales with the metal- 
licity of the outer stellar layers as (Z/Z©)" with a be- 
tween 0.5 and 0.8 (Kudritzki et al. [TM71 Vink & al. 
2001). Thus lowering the metaUicity by a factor 200 000 
(as would be the case for obtaining the metallicity of the 
Christlieb star) would thus lower the mass loss rates by 
a factor 450, or even by a greater factor if the metal- 
licity dependence becomes stronger at lower Z, as sug- 
gested by Kudritzki H2002|l . Now since metal-poor stars 
lose little mass, they also lose little angular momentum (if 
rotating), so they have a greater chance of reaching the 
break-up limit during the Main Sequence phase (see for 
instance Fig. 9 in Meynet & Maeder 2003)- At break-up, 
the outer stellar layers become unbound and are ejected 
whatever their metallicity. The break-up is reached more 
easily when we take into account that massive rotating 
stars have polar winds as shown by Owocki et al. H199(j|l 
and Maeder H1999|l . When most of the mass is lost along 
the rotational axis, little angular momentum is lost. 

Another way for rotation to trigger enhancements of 
the mass loss comes from the mixing induced by rotation. 
In general, rotational mixing favours the evolution into the 
red supergiant stage (see Maeder & Meynet EI^, where 
mass loss is higher. It also enhances the metallicity of the 
surface of the star and, in this way, boosts the radiatively 
driven stellar winds (see below). 

Thus there are very good reasons for exploring the ef- 
fects of rotation at very low metallicity, which we have 
attempted in this work. This was also the aim of the re- 
cent work by Marigo et al. H20()3|l . who compute Pop III 
massive stellar models with rotation, assuming solid-body 
rotation. In this very interesting piece of work, they mainly 
study the effects of reaching the break-up limit. However, 
since they did not include the rotational mixing of the 
chemical elements, they could not explore the effects of 
rotation on the internal chemical composition of the stars. 
Also, solid body rotation is just the extreme case of cou- 
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pling the internal rotation, which ignores the physics and 
timescales of the internal transport. In the present mod- 
els, the transport of both the angular momentum and the 
chemical species are treated in a consistent way, and, as 
we shall see, rotational mixing has very important conse- 
quences on both the stellar yields and the mass loss rates. 

In Sect. 2, we briefly recall the main physical ingredi- 
ents of the models. The evolutions of fast-rotating massive 
star models at very low Z are described in Sect. 3. The 
evolution of the internal chemical composition is the sub- 
ject of Sect. 4, while the ejected masses of various isotopes 
are presented in Sect. 5. The case of the CEMP stars is 
discussed in Sect. 6. Section 7 summarises the main results 
and raises a few questions to be explored in future works. 

2. Physical ingredients 

The computation of our models was done with the 
Geneva evolution code. The opacities were taken from 
Iglesias & Roger (1996) and complemented at low 
temperatures by the molecular opacities of Alexander 
( http : / / web . physics . twsu . edu/alex/wwwdra . htm). 
The nuclear reaction rates were based on the NACRE 
data basis ( Angulo & al. I1999|l . The treatment of rotation 
included the hydrostatic effects described in Meynet & 
Maeder H1997|l and the effects of rotation on mass loss 
rates according to Maeder & Meynet H2000(l . In partic- 
ular, we accounted for the wind anisotropies induced by 
rotation as in Maeder l|1999|l . Meridional circulation was 
implemented according to Maeder & Zahn H1998|l . but 
including the new Dh coefficient as described in Maeder 
H2003|l . Roughly compared to the old i?h, the new one 
tends to reduce the size of the convective core and to 
allow larger enrichment of the surface in CNO-processed 
elements. The reader is referred to these papers for a 
detailed description of the effects. The convective insta- 
bility was treated according to the Schwarzschild criterion 
without overshooting. The radiative mass loss rates are 
from Kudritzki & Puis when logTeff > 3.95 and 

from de Jager et al. ( 1988| otherwise. The mass loss rates 
depend on metallicity as M ~ (Z/Z©)^-^, where Z is 
the mass fraction of heavy elements at the surface of the 
star. As we shall see, this quantity may change during 
the evolution of the star. 

A specific treatment for mass loss was applied at break- 
up. At break-up, the mass loss rate adjusts itself in such a 
way that an equilibrium is reached between the two follow- 
ing opposite effects. 1) The radial inflation due to evolu- 
tion, combined with the growth of the surface velocity due 
to the internal coupling by meridional circulation, brings 
the star to break-up, and thus some amount of mass at 
the surface is no longer bound to the star. 2) By removing 
the most external layers, mass loss brings the stellar sur- 
face down to a level in the star that is no longer critical. 
Thus, at break-up, we should adapt the mass loss rates, 
in order to maintain the surface layers at the break-up 
limit. In practice, however, since the critical limit contains 
mathematical singularities, we considered that during the 



break-up phase, the mass loss rates should be such that the 
model stays near a constant fraction (for example, 0.98) 
of the limit. At the end of the MS phase, the stellar radius 
inflates so rapidly that meridional circulation is unable to 
continue to ensure the internal coupling, and the break-up 
phase ceases naturally. 

In this first exploratory work, we focused our atten- 
tion on stars with initial masses of 60 M© and 7 M© in 
order to gain insight into the properties of both massive 
and AGB stars at low Z . The evolution was computed un- 
til the end of the core helium burning phase (core carbon 
burning phase in the case of the 60 M© rotating model at 
Z — 10^^). Two metallicities were considered for the 60 
Mq models: Z = 10~^ and Z = 10"'^ Only this last metal- 
licity was considered for the 7 M© model. Of course we do 
not know if stars with Z — 10^^ have ever formed; how- 
ever, it is not possible at the present time to exclude such 
a possibility. Indeed it might be that the first star gener- 
ations produce very little amounts of heavy elements, due 
to the strong fallback of ejected material onto black holes 
at the end of their lifetimes. Moreover, as we shall see, the 
behaviour of the Z — 1Q~^ and 10~® massive star models 
are qualitatively similar, indicating that the evolutionary 
scenarios explored here might apply to a broad range of 
initial metallicities. 

The initial mixture of heavy elements was taken as 
equal to the one used to compute the opacity tables 
(Iglesias & Roger 119961 Weiss alpha-enhanced elements 
mixture). The initial composition for models at Z = 10~® 
is given in Table ^ The models at Z = 10~^ have the 
same initial mixture of heavy elements. More precisely, 
the mass fractions for all the isotopes heavier than ''He 
were multiplied by 10^ (=10-^10"^). 

Nothing is known on the rotational velocities of such 
stars. However, there are some indirect indications that 
stars at lower Z could have higher initial rotational veloc- 
ities: 1) realistic simulations of the formation of the first 
stars in the Universe show that the problem of the dissi- 
pation of the angular momentum is more severe at very 
low Z than at the solar Z. Thus these stars might begin 
their evolution with a higher amount of angular momen- 
tum (Abel et al. I2002|l . 2) There are some observational 
hints that the distribution of initial rotation might con- 
tain more fast rotators at lower Z (Maeder et al. I1999|l . 
3) Even if stars begin their life on the ZAMS with the 
same total amount of angular momentum for all metallic- 
ities, then the stars at lower metallicity rotate faster as a 
consequence of their smaller radii. 

The three arguments listed above would favour the 
choice of a higher value for the initial rotational veloc- 
ity than those adopted for solar models. To choose this 
value we proceeded in the following way. First we sup- 
posed that the stars begin their evolution on the ZAMS 
with approximately the same angular momentum content, 
whatever the metallicity. At solar metallicity, observation 
provides values for the mean observed rotational veloc- 
ity on the MS phase (around 200 km s~^ for OB stars). 
Stellar models allowed us to estimate the initial angular 
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Table 1. Initial abundances in mass fraction for models 
at Z = 10-^. 



Element 


Initial abundance 


H 


0.75999996 




0.00002554 


*He 


0.23997448 


12c 


7.5500C-10 


13q 


O.lOOOc-10 




2.3358C-10 




0.0092C-10 


16o 


67.100C-10 




0.0300C-10 


18(3 


0.1500C-10 


19p 


0.0020C-10 


^°Ne 


7.8368C-10 


21 Ne 


0.0200C-10 


^^Ne 


0.6306C-10 




0.0882C-10 


=*Mg 


3.2474C-10 


2^Mg 


0.4268C-10 


=«Mg 


0.4897C-10 


"Al 


0.1400C-10 




3.2800C-10 


56 Fo 


3.1675C-10 



momentum required to achieve such values (around 2.2- 
2.5 10^'^ g cm^ s~^). Adopting such an initial value of the 
angular momentum, we found that a 60 Mq stellar model 
at Z = 10~^ has a velocity on the ZAMS of 800 km s^^. 
This is the value of the initial velocity we adopt in the 
present work. 

3. Evolution of a massive rotating star at very low 
metallicity 

3.1. Rotation and mass loss during the main sequence 
phase 

Figure Q shows the evolutionary tracks during the main 
sequence (MS) phase for the 60 M0 stellar models at 
Z = 10~®. Table [5] presents some properties of the models 
at the end of the core H- and He-burning phases. From 
Fig. n we see that rotation produces a small shift of the 
tracks toward lower luminosities and TcS- This effect is 
due to both atmospheric distortions (note that surface- 
averaged effective temperatures are plotted in Fig. as 
explained in Meynet & Maeder I2000|l and to the low- 
ering of the effective gravity (see e.g. Kippenhahn and 
Thomas EZni Maeder and Pevtremann ll97()l Collins and 
Sonneborn ll977|l . The MS Hfetimc of the rotating model 
is enhanced by 11%. 

These results show that rotation does not affect the UV 
outputs of very metal-poor massive stars much (the UV 
outputs of the first massive star generations might con- 
tribute a lot to the reionization of the early Universe, see 
e.g. Madau 120035. Only if a significant fraction of primor- 
dial stars would rotate so fast that they follow the path 
of homogeneous evolution fMaeder Il987|l . could rotation 



increase the ionizing power. In that case, the star would 
remain in the blue part of the HR diagram and would have 
a much longer lifetime. 

Figure [3 shows the evolution of the ratio il/ilcrit at 
the surface during the MS phase. At Z ~ 10^^, the model 
with fini = 800 km reaches the break-up limit when 
the mass fraction of hydrogen at the centre — 0.40. 
The star stays at break-up for the remaining of its MS 
life with an enhanced mass loss rate. As a consequence, 
the model ends its MS life with 57.6 M0, having lost 4% 
of its initial mass. Despite the star staying in the vicinity 
of the break-up limit during an important part of its MS 
lifetime, it does not lose very large amounts of mass, due 
to the fact that only the outermost layers of the stars are 
above the break-up limit and are ejected. These layers 
have low density and thus contain little mass. 

A model with the same initial velocity, but with 
a metallicity three orders of magnitude higher, reaches 
the break-up limit very early in the MS phase, when 
Xc — 0.56. This comes from the fact that when the metal- 
licity increases, a given value of the initial velocity corre- 
sponds to a higher initial value of the Wini/^-'crit ratio. The 
model at Z = 10"^ ends its MS life with 53.8 Mq, having 
lost 10% of its initial mass. 

In order to discuss the effects of a change of rotation, 
it is interesting to compare this last result with the one 
obtained in Meynet & Maeder H2002|l for a 60 Mq model at 
Z — 10~^ with Vini= 300 km s~^. This last model reaches 
the break-up velocity much later, only when Xc — 0.01. 
At Z = 10~^, the velocity 300 km s~^ thus appears as 
the lower limit for the initial rotation, allowing a 60 M© 
star to reach the break-up limit during its MS phase. The 
60 Mq star ends its MS life with 59.7 Mq, having lost 
only 0.5% of its initial mass. Note that the 2002 models 
were computed with slightly different physical ingredients 
than those used to compute the models discussed in this 
paper (different expression for Dh and for the mass loss 
rates). However, at these very low metallicities, radiatively 
driven winds remain quite modest, and the transport of 
the angular momentum, mainly driven by the meridional 
circulation, does not depend much on the expression of 

During the MS phase, the surface of the rotating stars 
is enriched in nitrogen and depleted in carbon as a result 
of rotational mixing. Figure 13 shows that the N/C ratios 
are enhanced by more than two orders of magnitude at 
the end of the H-burning phase. More precisely, at the 
surface of the Z = 10~^ model, nitrogen is enhanced by a 
factor 27 and carbon decreased by a factor 12. All other 
physical ingredients being the same, one also sees that 
the model with the lowest metallicity is also the one with 
the greatest surface enrichments. This well agrees with 
the trend already found in our previous works (Maeder 
& Meynet EHHl Meynet & Maeder EEEj), which results 
from the steep gradients of angular velocity that build up 
in very metal-poor stars that favours shear mixing. Let us 
emphasise, however, that during the MS phase, the arrival 
of CNO-processed material at the surface does not change 
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Table 2. Properties of the stellar models at the end of the H- and He-burning phases. Afi^i is the initial mass, Z the 
initial metallicity, Ujni the initial velocity, v the mean equatorial rotational velocity during the MS phase defined as in 
Meynet & Maeder (j2()()()|l . in the H-burning lifetimes, M the actual mass of the star, v the actual rotational velocity 
at the stage considered, the helium surface abundance in mass fraction. N/C and N/0 are the ratios of nitrogen 
to carbon, respectively, of nitrogen to oxygen at the surface of stars in mass fraction; C, N, O are the abundances of 
carbon, nitrogen, and oxygen at the surface in mass fractions. The numbers in parentheses indicate the power of ten, 
i.e., 7.54(-10)=7.54 xlO"!". 













End of H- 


-burning 












End of He 


—burning 








Z 


^ini 


V 


tH 


M 


V 




N/C 


N/O 


tHc 


M 


V 




C 


N 


O 


Ma 




km 


km 


Myr 


Mq 


km 








Myr 


Mq 


km 










60 


10"** 








3.605 


59.817 





0.24 


0.31 


0.03 


0.292 


59.726 





0.24 


7.54(-10) 


2.34(-10) 


6.71(-9) 


60 


10"** 


800 


719 


4.004 


57.624 


591 


0.27 


103 


9.77 


0.522 


23.988 


0.02 


0.76 


1.97(-4) 


1.02(-2) 


2.85(-4) 


60 


10"^ 


800 


636 


4.441 


53.846 


567 


0.34 


40 


0.82 


0.544 


37.280 


0.57 


0.80 


5.06(-5) 


2.07(-3) 


4.24(-5) 




4.9 4.85 4.8 4.75 4.7 

log T-,,, [K] 

Fig. 1. Evolutionary tracks in the HR diagram for 60 Mq 
stellar models a.t Z = 10^*. 



the metallicity of the outer layers. Indeed, rotational mix- 
ing brings nitrogen to the surface but depletes carbon and 
oxygen, keeping the sum of CNO elements constant, and 
the metallicity as well. Thus during the MS phase, the en- 
hanced mass loss rates undergone by the rotating models 
are entirely due to the mechanical effect of the centrifugal 
force. As we shall see, this is no longer the case during the 
core He-burning phase. 

3.2. Rotation and mass loss during the post 
main-sequence pliases 

Rotational mixing that occurs during the MS phase and 
still continues to be active during the core He-burning 
phase deeply modifies the internal chemical composition 
of the star. This has important consequences during the 
post-MS phases and deeply changes the evolution of the 
rotating models with respect to their non-rotating coun- 
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Fig. 2. Evolution of fi/ficrit at the surface of 60 M© mod- 
els ai Z = 10^* (continuous line) and Z = 10^^ (upper 
dashed line) with z;ini= 800 km s~^. The case of the 60 M© 
model aX Z = 10~^ with Vin\= 300 km s~^ from Meynet 
& Maeder H2002|l is also shown (lower dashed line). 

terparts. Among the most striking differences, one notes 
the following: 

— 1) Rotation favours redwards evolution in the HR di- 
agram as was already shown by Maeder & Meynet 
©Ol), and as illustrated in Fig. ^ One sees that the 
non-rotating model remains on the blue side during 
the whole core He-burning phase, while the 800 km 
s^^ model aX Z — 10~^ starts its journey toward the 
red side of the HR diagram early in the core helium 
burning stage, when Yc ~ 0.67 {Y^ is the mass frac- 
tion of helium at the centre of the star model). The 
same is true for the corresponding model a,t Z = 10~^. 
Let us recall that this behaviour is linked to the rapid 
disappearance of the intermediate convective zone as- 
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Fig. 3. Evolution as a function of log Tea of the excess at 
the surface for the ratio N/C (expressed in dex) compared 
to the initial ratio. N and C are the abundances of nitrogen 
and carbon at the surface of the star. 
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Fig. 4. Evolution of logTcff as a function of Vc, the mass 
fraction of ^He at the centre, for a non-rotating (dashed 
line) and rotating (continuous line) 60 model at 
Z= 10-8. 



sociated to the H-burning shell (see Fig. [S] and Maeder 
& Meynet EnnHl- 

2) Redwards evolution enhances the mass loss. In 
the cases of our 60 M© stellar models, it brings the 
stars near the Humphreys-Davidson limit, i.e., near 
logL/L© — 6 and Tcff in a broad range around 10^ 
K. Near this limit, the mass loss rates (here from 
de Jager et al. 11988(1 become very important. For in- 
stance, the model represented in the left panel of Fig.|Sl 
(logL/L© = 6.129, logTcff = 4.243) is stih far to the 
left hand side of the Humphreys-Davidson limit. Its 
mass loss rate is log(— M) = —5.467, where M is ex- 
pressed in Mq per year. The model in the right panel 
(logL/L0 = 6.145, logTcff = 3.853) is in the vicinity 
of the Humphreys-Davidson limit. Its mass loss rate 



is equal to -4.616, i.e., more than seven times higher. 
During this transition, the overall metallicity at the 
surface does not change and remains equal to the ini- 
tial one (here Zi„i = 0.00001). We observe a similar 
transition in the case of the Z — 10"^ stellar model. 

— 3) During the core He-burning phase, primary nitro- 
gen is synthesized in the H-burning shell, due to the 
rotational diffusion of carbon and oxygen produced in 
the helium core into the H-burning shell (Meynet & 
Maeder 2002). This is illustrated well in Fig. [H for the 
Z = 10^^ rotating model and in Fig.|Blfor the model 
at Z ^ 10-*. 

— 4) In contrast to what happens during the MS phase, 
rotational mixing during the core He-burning phase 
induces large changes in the surface metallicity. These 
changes occur only at the end of the core He-burning 
phase, although the conditions for their apparition re- 
sult from the mixing that occurs during the whole core 
He-burning phase. Indeed, rotational mixing progres- 
sively enriches the outer radiative zone in CNO ele- 
ments, thus enhancing its opacity slowly. When, in the 
Z = 10^8 stellar model, the abundance of nitrogen in 
the outer layers becomes approximately 10^* in mass 
fraction (i.e., has increased by two orders of magni- 
tude with respect to the initial value), these outer lay- 
ers become convective. The outer convective zone then 
rapidly deepens in mass and dredges up newly synthe- 
sized elements to the surface. From this stage onwards, 
the surface metallicity increases in a spectacular way, 
as can be seen in Fig. El For instance, the rotating 60 
Mq aX Z — 10^8 has a surface metallicity of 10^^ at 
the end of its lifetime, i.e., similar to that of the Large 
Magellanic Cloud ! 

— 5) The consequence of such large surface enrichments 
on the mass loss rates remains to be studied in de- 
tail using models of stellar winds with the appropriate 
physical characteristics (position in the HR diagram 
and chemical composition). In the absence of such so- 
phisticated models, we applied the usual rule here, 
namely M{Z) {Z/ZqYI'^M{Zq), where Z is the 
metallicity of the outer layers. With this prescription, 
the surface enhancement of the metallicity is responsi- 
ble for the large decrease in the stellar mass that can 
be seen in Fig. El 

— 6) During the late stages of the core helium-burning 
phase, as a result of mass loss and mixing, the star 
may evolve along a blue loop in the HR diagram (see 
Fig. 01). When the star evolves bluewards, the global 
stellar contraction brings the outer convective zone, 
which evolves like a solid body rotating shell to break- 
up (Heger & Langer [l998|l . At this stage of the evo- 
lution, the luminosity is not far from the Eddington 
limit and the star may reach the fiF-limit (Maeder & 
Meynet I2()()()|l . This multiplies the mass loss rates by 
very large factors. 

— 7) During the last 24 000 years of its lifetime, the 
model presents abundance patterns characteristic of 
WNL stars at its surface. 
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— 8) As a result of mixing and mass loss rates, the du- 
ration of the core He-burning phase is much longer in 
the rotating model. The present 60 Mq model with 
fini = 800 km s~^ at Z = 10~^ has a helium-burning 
lifetime that is '^80% longer than the corresponding 
lifetime of the non-rotating model (see Sect. 4 below 
for more explanations). 

The different effects described above are all due to the 
mixing induced by rotation and they all tend to enhance 
the quantity of mass lost by stellar winds. The rotating 
60 M0 at Z = 10^^ loses about 36 M0 during its lifetime. 
About 2 Mq arc lost due to break-up during the MS phase, 
~ 3 Mq are lost when the star is in the red part of the HR 
diagram (with surface metallicity equal to the initial one), 
27 Mq are lost due to the effect of the enhancement of the 
surface metallicity, the remaining 4 Mq are lost when the 
star evolves along the blue loop and reaches the riF-limit. 
One sees that, by far, the most important effect is due to 
the increase in the surface metallicity. 

Paradoxically the corresponding model at higher 
metallicity {Z = 10~^) loses less mass (a little less than 
40% of the total mass). This can be understood from the 
following facts: first less primary nitrogen is synthesized 
due to slightly less efficient chemical mixing when the 
metallicity increases, thus the surface does not become 
as metal rich as in the model at Z ^ 10^^. Second and for 
the same reason as above, the outer convective zone does 
not deepen as far as in the more metal-poor model. These 
two factors imply that the maximum surface metallicity 
reached in this model, which is equal to 0.0025, is about 
a factor 4 below the one reached by the Z = 10^^ model. 
Finally, the blue loop does not extend that far into the blue 
side, and the surface velocity always remains well below 
the break-up limit during the whole blueward excursion. 

In order to investigate to what extent the behaviour 
described above depends on the physical ingredients of 
the model, we compare the present results with those of 
a rotating model (fini= 800 km s~^) of a 60 Mq star 
at Z = 10~^ with a different prescription for the mass 
loss rates (Vink et al. I MM instead of Kudritzki 

& Puis I20dn|l . with the Ledoux criterion instead of the 
Schwarzschild one for determining the size of the con- 
vective core, with a core overshoot of a = 0.2 Hp and 
the old prescription for the horizontal diffusion coefficient 
Dh- This model is described in Meynet et al. (2005). In 
this case, the outer convective zone deepens farther into 
the stellar interior and thus produces a greater enhance- 
ment of the surface metallicity (the same order as the 
one we obtained in the present Z = 10~* 60 Mq model). 
Higher enhancements of the surface metallicity then in- 
duces greater mass loss by stellar winds. More important 
than these differences, however, we shall retain here that 
the results are qualitatively similar to those obtained in 
our previous models. In particular, the mechanism of sur- 
face metallicity enhancement occurs in both models and 
appears to be a robust process. 



3.3. Do very metal-poor, very massive stars end their 
lives as pair-instability supernovae 7 

Might the important mass loss undergone by the rotat- 
ing models prevent the most massive stars from going 
through pair instability ? According to Heger & Woosley 
H2002|l , progenitors of pair-instability supernovae have he- 
lium core masses between '^64 and 133 Mq. This corre- 
sponds to initial masses between about 140 and 260 Mq. 
Thus the question is whether stars with initial masses 
above 140 Mq can lose a sufficient amount of mass to 
have a helium core that is less than about 64 Mq at the 
end of the core He-burning phase. From the values quoted 
above, it would imply the loss of more than (140-64)=76 
Mq, which represents about 54% of the initial stellar mass. 
From the computation performed here, one can expect 
that such a scenario is possible, where a 60 Mq loses more 
than 60% of its initial mass. However, more extensive com- 
putations are needed to check whether the rotational mass 
loss could indeed prevent the most massive stars from go- 
ing through this pair instability. Were this the case, it 
would explain why the nucleosynthetic signature of pair- 
instability supernovae is not observed in the abundance 
pattern of the most metal-poor halo stars observed up 
to now. At least this mechanism could restrain the mass 
range for the progenitors of pair-instability supernovae, 
pushing the minimum initial mass needed for such a sce- 
nario to occur to higher values. Let us also note that the 
luminosity of the star comes nearer to the Eddington limit 
when the initial mass increases. When rotating, such stars 
will then encounter the ilF- limit (Maeder & Mevnet F2000|l 
and very likely undergo strong mass losses. 

4. Evolution of the interior chemical composition 

As discussed above, rotational mixing changes the chem- 
ical composition of stellar interiors in an important way. 
This is illustrated well by Figs. El and which show the 
internal chemical composition of our rotating and non- 
rotating 60 Mq stellar models at four different stages 
at the end of the core He-burning phase (models at 
Z = 10-8). 

Comparing panels a of Figs. Eland [7| one sees that a 
large convective shell is associated to the H-burning shell 
in the rotating model, while such a shell is absent in the 
non-rotating model. This contrasts with what happens at 
the beginning of the core He-burning phase, where the in- 
termediate convective zone associated to the H-burning 
shell was absent in the rotating model (or at least much 
smaller), while in the non-rotating model, the intermedi- 
ate convective zone was well-developed (see above). Why 
is there this difference between the beginning and the end 
of the core He-burning phase ? At the beginning of the 
core He-burning phase, the disappearance of the interme- 
diate convective shell was a consequence of the rotational 
mixing that operated during the core H-burning phase and 
that brought some freshly synthesized helium into this re- 
gion. More helium in this region means less hydrogen and 
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Fig. 5. Chemical composition of a 60 M© stellar model at Z = 10~^ with Wini= 800 km s~^ when it evolves from the 
blue to the red part of the HR diagram. The model shown in the left panel has logL/L© = 6.129 and logTeff = 4.243; 
in the middle panel, it has logL/L© — 6.130 and logToff = 4.047; in the right panel, it has logL/L© = 6.145 and 
logTeff = 3.853. 



also some decrease in the opacity, both of which inhibit the 
development of convection (cf . Maeder & Meynet 12001(1 . 
Now, at the end of the core He-burning phase, we have 
He-burning products that are brought into the H-burning 
shell. These products, mainly carbon and oxygen, act as 
catalysts for the CNO cycle and make the H-burning shell 
more active, thus favouring convection. This mechanism 
enriches this zone not only in primary ^''N, but also in 
primary "'^^C. Looking at the H-rich envelope in the non- 
rotating model at the same stage, one sees that all these 
elements have much lower abundances. Actually they fall 
well below the minimum ordinate of the figure. 

If one now compares the chemical composition of the 
CO cores when Yc ^ 0.11 (see panels a of Figs. El and [T)), 
one notes the following points. First, the abundances in 
^^C, '^^O and ^"Ne are approximately equal in both the ro- 
tating and non-rotating models. This comes from the fact 
that the CO core masses are approximately the same in 
both models. On the other hand in the rotating model, the 
abundance of ^^Ne is greatly enhanced, as the abundances 



of ^^Mg and ^^Mg. The abundance of ^^Ne results from 
the conversion of primary ^"^N, which has diffused into 
the He-burning core. Thus the resulting high abundance 
of ^^Ne is also of primary origin. The isotopes of magne- 
sium are produced by the reactions ^^Ne(a, 7)^^Mg and 
^^Ne(Q!, n)^^Mg. Their high abundances also result from 
primary nitrogen diffusion into the He-core. 

The CO-core mass (cf. Table IS)) in the rotating model 
is slightly smaller than in the non-rotating one. This con- 
trasts to what happens at higher metallicity, where rota- 
tion tends to increase the CO-core mass (see Hirschi et al. 
W^. A gain, this results from the mechanism of primary 
nitrogen production, which induces a large convective zone 
associated to the H-burning shell, which then prevents this 
shell from migrating outwards and, thus, the CO core from 
growing in mass. Let us recall that in rotating models at 
solar metallicity, there is no primary nitrogen production 
due to the less efficient mixing at higher metallicity (see 
Meynet & Maeder [2002(1 : thus there is no increase in the 
H-burning shell activity. 
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Fig. 6. Variations of the abundances (in mass fraction) as a function of the Lagrangian mass within a 60 Mq star 
with Vhii— 800 km s~^ and Z — 10~*. The four panels show the chemical composition at four different stages at the 
end of the core He-burning phase: in panel a) the model has a mass fraction of helium at the centre, Yc= 0.11 and an 
actual mass M= 54.8 Mq - b) Yc= 0.06, M= 48.3 Mq - c) Yc= 0.04, M= 31.5 Mq - d) End of the core C-burning 
phase, M= 23.8 Mq. The actual surface metallicity Zsmf is indicated in each panel. 



In panel b of Fig. El as explained in the previous sec- 
tion, one sees the outer convective zone extending inwards 
and bringing CNO elements to the surface. In panels c and 
d, mass-loss efficiently removes these outer layers. At the 
corresponding stages in the non-rotating model, the outer 
envelope is not enriched in heavy elements and keeps its 
mass. 

At the end of the He-burning phase (see panels d), 
the abundance of ^^C is significantly smaller in the ro- 
tating model than in the non-rotating one. At the same 
time, the abundances of ^°Ne and ^''Mg are significantly 
greater. This is a consequence of helium diffusion into the 
He-core at the end of the He-burning phase. Let us recall 
that ^^C is destroyed by alpha capture (to produce ^®0), 
while ^"Ne and ^''Mg are produced by alpha captures on, 
respectively, ^^O and ^°Ne. For what concerns the other 
isotopes of neon and magnesium, one sees that in the ro- 
tating models, much higher abundances of ^^Mg and ^^Mg 
are reached due to the transformation of the ^^Ne at the 



end of the core helium-burning phase. The neutrons lib- 
erated by the ^■^Ne(Q;,n)^^Mg reaction can be captured by 
iron peak elements, producing some amount of s-process 
elements (see e.g Baraffe et al. I1992f) . In view of the im- 
portant changes to the interior chemical composition due 
to rotation, there is good chance that the s-process in the 
present rotating massive star models is quite different from 
the one obtained in non-rotating models. This will be ex- 
amined in later papers. 

5. Chemical composition of the winds and of the 
supernova ejecta 

5.1. Wind composition 

Let us first discuss the chemical composition of the winds. 
The total mass lost, as well as the quantities of various 
chemical elements ejected by stellar winds, are given in 
Tables 01 & H The models at Z = 10"^ were computed 
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Fig. 7. Same as Fig. Elfor a 60 Mq star with Vini= km 
s~^ and Z ~ 10~*. The four panels shows the chemical 
composition at four different stages at the end of the core 
He-burning phase: Panel a) Yc= 0.12, M= 59.74 M© - b) 
Fc= 0.06, M= 59.74 M© - c) Y^^ 0.03, Af= 59.73 M© 
- d) Yc= 0.0, M= 59.73 Mq. The surface metalHcity is 
equal to 10~^ at the four evolutionary stages. 

with an extended nuclear reaction network including the 
Ne-Na and Mg-Al chains, which is why in Table 0] the 
wind-ejected masses of these elements can be indicated. 
The stellar yields - i.e., the mass of an isotope newly syn- 
thesized and ejected by the star - can be obtained by sub- 
tracting the mass of the isotope initially present in that 
part of the star^ from the ejected masses given in Tables O 
&I1 

According to Table O the non-rotating model ejects 
only half a percent of its total mass through stellar winds, 
which is completely negligible. Moreover, this material 
has exactly the same chemical composition as does the 
protostellar cloud from which the star formed. If, at the 
end of its lifetime, all the stellar material is swallowed by 
the black-hole resulting from the star collapse, the nu- 
cleosynthetic contribution of such stars would be zero. 
In contrast, the rotating models lose more than 60% of 
their initial mass through stellar winds. This material is 
strongly enriched in CNO elements. Even if all the final 
stellar mass is engulfed into a black-hole at the end of the 
evolution, the nucleosynthetic contribution of such stars 
remains quite significant. As already noted above, the cor- 
responding model at Z = 10^^ loses less mass by stellar 
winds (see Table ^J. However, the amounts of mass lost 
remain large and they present strong enrichments in CNO 

^ Tliis quantity may be obtained by multiplying the initial 
abundance of the isotope considered (given in Table0 by rriej. 



elements, as in the case of the Z ~ 10 ^ rotating model. 
Also 2^Na and ^^Al are somewhat enhanced in the wind 
material. 

Other striking differences between the rotating and 
non-rotating models concern the ^^C/^'^C, N/C, and N/0 
ratios (see Tables |31 & 0J . The wind of the non- rotating 
model shows solar ratios {^^G/^^G=73, N/C=0.31, and 
N/O=0.03 in mass fractions). The wind of rotating mod- 
els is characterised by very low ^^C/^^C ratios, around 4 
- 5 (close to the equilibrium value of the CN cycle) and 
by very high N/C (between about 3 and 40) and N/0 ra- 
tios (between 1 and 36). Thus wind material presents the 
signature of heavily CNO-processed material. 

5.2. Total ejecta composition (wind and supernova 
ejecta ) 

In order to estimate the quantity of mass lost at the time 
of the supernova explosion (if any) , it is necessary to know 
the mass of the remnant. This quantity is estimated with 
the relation of Arnett 1)1991(1 between the mass of the rem- 
nant and the mass of the carbon-oxygen core. The masses 
of the different elements ejected are then simply obtained 
by integrating their abundance in the final model between 
TOrom (see Tables 01 & ^ and the surface. Since the evo- 
lution of the present models was stopped before the pre- 
supernova stage was reached, the masses of ^^C and ^^O 
obtained here might still be somewhat modified by the 
more advanced nuclear phases. 

How does the contribution of the two models (rotating 
and non-rotating) at Z — 10~* compare when both the 
wind and the supernova contribute to the ejection of the 
stellar material? First, one sees that the total mass ejected 
(through winds and supernova explosion) is very similar 
(on the order of 54-55 Mq), due to the fact that the two 
models have similar CO core masses. Second, one sees that 
the amount of '^He ejected by the rotating model is slightly 
higher, whereas the amount of ^^C is lower due to the ef- 
fect discussed above (a-captures at the end of the core 
He-burning phase). The quantity of ^^O ejected is similar 
in both models. Third, the most important differences be- 
tween the rotating and non-rotating models occur for ^^C, 
-'^'^N, ^^O, and ^^O. The abundances of these isotopes are 
increased by factors between 10*-10^ in the ejecta of the 
rotating model. The first three isotopes are produced in 
the H-burning shell (CNO cycle) and are mainly ejected by 
the winds, while the last one, produced at the interface be- 
tween the CO-core and the He-burning shell, is ejected at 
the time of the supernova explosion. Fourth, rotation also 
deeply affects the ratios of light elements in the ejected 
material (see Tables O & 0J) . The effects of rotation are 
qualitatively similar to those obtained when comparing 
the composition of the wind material of rotating and non- 
rotating stellar models. Rotation decreases the ^-^C/^^C 
ratios from 3.5 x 10® in the non-rotating case to 311 in 
the rotating case, while it increases the N/C and N/0 ra- 
tios, which have values of ~10"^ and 10^® respectively, 
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when Uini = km s ^, and of 0.5 and 0.02 when Vini = 
800 km s-i. 

In the ejecta of rotating models, composed of both 
wind and supernova material, the ^^C/^'^C ratio is higher 
than in pure wind material, and the N/C and N/O ratios 
are smaller. This comes from the fact that the supernova 
ejecta are rich in helium-burning products characterized 
by a very high ^^C/^^C ratio and by very low N/C and 
N/O ratios. 

At this point we can ask if the rotating star had lost no 
mass through stellar winds, and if all the stellar material 
were ejected at the time of the supernova explosion, would 
the composition of the ejecta be different with respect to 
the case discussed above, where part of the material is 
ejected by the winds and part by the supernova explo- 
sion. Let us recall that stellar winds remove layers from 
the stars at an earlier evolutionary stage than do super- 
nova explosions. If some of these layers, instead of being 
ejected by the winds, had remained locked inside the star, 
they would have been processed further by the nuclear 
reactions. Thus their composition at the end of the stel- 
lar evolution would be different from the one obtained if 
they had been ejected at an earlier time by the winds. 
Obviously for such differences to be important, mass loss 
must remove the layers at a sufficiently early time. If it 
does so only at the very end of the evolution, there would 
be no chance for the layers to be processed much by the 
nuclear reactions, and there would be no significant differ- 
ence whether the mass were ejected by the winds or by the 
supernova explosion. Actually, this is what happens in our 
rotating models. As indicated above, the mass is removed 
at the very end of the He-burning phase, and only material 
from the H-rich envelope is ejected. Thus, if this material 
were ejected only at the time of the supernova explosion, 
it would have kept the same chemical composition as the 
one in Table [31 

As a result, the chemical composition of the ejecta 
(wind and supernova) does not depend much on the mass 
loss, but is deeply affected by rotation. However, the stel- 
lar winds may of course be of primary importance if the 
whole final mass of the star is swallowed in a black-hole 
at the end of the evolution. In that case, the star will con- 
tribute to the interstellar enrichment only by its winds. 

Comparing the data given in the right part of Table 13 
with the left part of Table ^ one can see the effect 
produced by an increase in the initial metallicity by 
three orders of magnitude (all other things being equal). 
Interestingly, we see that the differences between the two 
models are in general much smaller than those between 
the rotating and non-rotating model at a given metallic- 
ity. The total ejected mass, and the masses of ''He, ^^C, 
"'^O are similar within factors between 0.8 and 1.4. The 
quantities of ^^N and ^^C are within an order of magni- 
tude, and the masses of ^''O and ^^O differ by a factor 3. 
We are thus far from the factors 10^-10^ between the re- 
sults of the rotating and non-rotating models at Z = 10~* 
! The effects of rotation at extremely low metallicity are 
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Table 3. Helium-, CO-core mass and mass of the rem- 
nants (respectively rua, nico, and mrom) of 60 M0 stellar 
models with and without rotation at Z = 10~*. The total 
mass ejected (woj) and the mass ejected of various chem- 
ical species {■m{Xi)) are given in solar masses. The values 
of some isotope ratios (in mass fractions) are also indi- 
cated. The case of matter ejected by stellar winds only is 
distinguished from the case of matter ejected by both the 
stellar winds and the supernova explosion. 





Mini/Mg Z Uini 


Mini/Mg Z Ui„i 




60 10"** 


60 IQ-** 800 


mco 


23.08 
21.61 
6.65 


23.83 
18.04 
5.56 




Mass ejected 


Mass ejected 




WIND SN+WIND 


WIND SN+WIND 


"loj 

mC*Ho) 
mC^C) 
m("C) 
m(i*N) 

mCO) 
mCO) 


0.28 53.35 
6.62e-02 19.58 
2.08e-10 2.066 
2.84e-12 5.84e-09 
6.44e-ll 1.90e-07 
1.85e-09 12.61 
8.27e-13 6.39e-10 
4.14e-12 6.53e-10 


36.17 54.44 
21.46 23.85 
4.78e-03 4.26e-01 
1.25e-03 1.37e-03 
1.97e-01 2.20e-01 
6.08e-03 13.54 
7.50e-06 8.60e-06 
1.58e-08 4.44e-03 




Isotopic ratios 


Isotopic ratios 


N/C 
N/O 


73.24 3.54e+08 
0.31 9.20e-08 
0.03 1.51e-08 


3.82 311 
41.2 0.52 
32.4 0.02 



much larger than the effects of a change in the initial Z 
content. 

The results given on the right side of Table 0] corre- 
sponds to the model described in Meynet et al. (,.2005j) . It 
differs from the present models by the mass loss and mix- 
ing prescription (see Sect. 3.2). As already emphasized 
above, the results are qualitatively very similar. However, 
quantitatively, they present some differences. For instance, 
the quantity of ^^C in the model presented on the right 
side of Table ^ is larger by a factor 4 compared to the 
value given on the left side of the same Table. The right 
model presents a smaller helium core, an effect mainly 
due to higher mass loss rates. This favours larger ejec- 
tions of carbon by the winds and also by the supernova, 
since smaller helium cores lead to higher C/0 ratios at 
the end of the helium-burning phase. In the right model, 
the quantity of ^^O is decreased by about 4%. The ejected 
masses of ^^C and ^"^O are increased by factors of 176 and 
55, respectively. The masses of the other isotopes differ by 
less than an order of magnitude. 



12 



G. Meynet, S. Ekstrom, & A. Maeder: The early star generations 



Table 4. Same as Table 13 for rotating stellar models at 
Z — 10"^. The two 60 Mq models were computed with 
different physical ingredients, see text. 



m(*Hc) 
m(i=C) 
m("C) 
m("N) 
m("N) 
m(i<'0) 
m("0) 
m(i«0) 
m(i^F) 
m(=*''Ne) 

m(^^Na) 
m('*Mg) 
m('"Mg) 
m(^<'Mg) 
m(2^Al) 



N/C 
N/O 



60 10"-' 800 



36.90 
28.60 
8.69 



Mass ejected 
WIND SN+WIND 



23.10 
12.20 
3.34e- 
6.81e- 
9.78e- 
3.21e- 
2.72e- 
4.59e- 
2.82e- 
1.95e- 
7.64e- 
1.84e- 
4.80e- 
1.22e- 
4.41e- 
3.95e- 
5.68e- 
5.49e- 



04 
05 
03 
07 
04 
07 
08 
09 
06 
08 
07 
06 
06 
07 
07 
06 



51.31 

18.99 

5.84O-01 

1.47e-04 

2.51e-02 

2.15e-06 

18.12 

2.40e-06 

1.34e-03 



4.90 
29.3 
36.0 



Isotopic ratios 

3970 
0.04 
0.001 



Mi„i/Mo 



60 



10" 



800 



30.69 
27.95 
8.50 



Mass ejeeted 
WIND SN+WIND 



29.31 
12.60 
1.45e- 
3.81e- 
4.29e- 
1.55e- 
3.29e- 
2.78e- 
1.63e- 
l.lOe- 
1.29e- 
6.99e- 
3.35e- 
5.61e- 
6.21e- 
6.96e- 
3.07e- 
7.75e- 



02 
03 
02 
06 
02 
05 
08 
08 
05 
08 
05 
06 
06 
07 
06 
06 



51.50 
14.58 
2.46 

2.58e-02 

1.87e-01 

1.68e-05 

17.32 

1.32e-04 

2.07e-04 



3.81 
2.96 
1.30 



Isotopic ratios 

95.3 
0.08 
0.01 



6. Link with the extremely metal-poor C-rich 
stars 

6.1. Observations and existing interpretations 

Spectroscopic surveys of very metal-poor stars (Beers et 
al. Beers [T^M Christlieb |2)03j have shown that 

CEMP stars account for up to about 25% of stars with 
metallicities lower than [Fe/H]'^ —2.5 (see e.g. Lucatello 
et al. EOnil)- A star is said to be C-rich if [C/Fe]> 1. A 
large proportion of these CEMP stars also present en- 
hancements in their neutron capture elements (mainly s- 
process elements). A few of them also appear to exhibit 
large enhancements in N and O. The most iron-deficient 
stars observed so far are CEMP stars. These stars are HE 
0107-5240, a giant halo star, and HE 1327-2326, a dwarf 
or subgiant halo star. The star HE 0107-5240 ([Fe/H]=- 
5.3) presents the following CNO surface abundances: 
[C /Fe]= 4.0, [N/Fe]=2 .3, an d [0/Fe]=2.4 (Christlieb et 
al. EODl Bessell et ah OHIMjl . The ratio ^^C/^^C has also 
been tentatively estimated by Christlieb et al. H2004[l . who 



suggest a value of about 60, but with a great uncertainty. 
They can, however, rule out a value inferior to 50 (let us re- 
call that the solar ratio is ~73). The star HE 1327-2326 has 
[Fe/H]=-5.4 and CNO surface abundances: [C/Fe]=4.1, 
[N/Fe]=4.5, [0/Fe]< 4.0 (Frebel et aL lM75|l . 

The origin of the high carbon abundance is still an 
open question and various scenarios have been proposed: 

1. The primordial scenarios: in this case the abun- 
dances observed at the surface of CEMP stars are the 
abundances of the cloud from which the star formed. 
The protostellar cloud was enriched in carbon by one 
or a few stars from a previous generation. For instance, 
Umeda and Nomoto (2003) propose that the cloud 
from which HE 0107-5240 formed was enriched by the 
ejecta of one Pop HI 25 M© star, which had exploded 
with low explosion energy (on the order of 3 x 10^° erg) 
and experienced strong mixing and fallback at the time 
of the supernova explosion. The mixing is necessary to 
create the observed high-level enrichments in light el- 
ements, and the fallback is necessary to retain a large 
part of the iron peak elements. Limongi et al H2()()3|l 
suggest that the cloud was enriched by the ejecta of 
two supernovae from progenitors with masses of about 
15 and 35 Mq. 

2. The accretion/ mass transfer scenarios: some au- 
thors have proposed that this particular abundance 
pattern results from accretion of interstellar material 
and from a companion (for instance an AGB star, as 
proposed by Suda et al. i20041 . As far as the nucleosyn- 
thetic origin is concerned, this scenario is not funda- 
mentally different from the first one. 

3. The in situ scenarios: finally, some authors have 
explored the possibility that the star itself could have 
produced the particular abundance pattern seen at its 
surface (Picardi et al. i2004j) . The overabundance of ni- 
trogen might easily be explained in the frame of this 
scenario, if the star had begun its evolution with the 
high carbon and oxygen overabundance. In fact we did 
perform a test calculation of a non-rotating 0.8 M© 
stellar model at [Fe/H]=-5.3 with an initial value of 
[C/Fe] and [0/Fe] equal to 4.0 and 2.3, respectively, 
i. e. equal to the abundances observed at the surface of 
HE 0107-5240. We found that, when the star reaches 
the value of the effective temperature (Tes = 5100± 
150 K) and of gravity (log 5 = 2.2 ± 0.3) of HE 0107- 
5240 (Christlieb et al. 'SUOT), the surface nitrogen en- 
richment is well within the range of the observed val- 
ues. However, it appears difficult to invoke similar pro- 
cesses to explain the high carbon and oxygen enhance- 
ments (see Picardi et al. EE 



6.2. No "in situ" CN production 

An abundance pattern typical of CEMP stars has been 
observed at the surface of non-evolved stars (Norris et 
aL lTWI Plez & Cohen Frebel et al.[lD05,). Among 

the most recent observations, let us mention the subgiant 
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or dwarf star HE 1327-2326 (Frebel et al. lMTSl see above) 
and the dwarf star G77-61 (Plez & CohenP05 ). The ini- 
tial mass of G77-61 is estimated to be between 0.3 and 0.5 
Mq, and it has an [Fc/H]=-4.03, [C/Fe]=2.6, [N/Fe]=2.6, 
and a ^^C/^^C ratio of 5±1. In this case, there is no way 
for the star, which burns its hydrogen through the pp 
chains, to produce nitrogen. There is even less possibility 
of producing surface enhancements of carbon and oxygen. 
Therefore, the "in situ" scenario can be excluded, at least 
for this star. In that case, only the first and second sce- 
narios are possible. The same is true for explaining the 
very high overabundances of carbon and nitrogen at the 
surface of HE 1327-2326. 

The values observed at the surface of non-evolved and 
evolved stars are shown in Fig. |S1 In the second case, the 
surface may have been depleted in carbon and oxygen 
and enriched in nitrogen due to the dredge-up that oc- 
curs along the red giant branch. On the other hand, in 
the case of the non-evolved stars, as explained above, this 
mechanism cannot be invoked, and the measured abun- 
dances reflect the abundances of the cloud that gave birth 
to the star. On the whole, the distributions of elements is 
similar for evolved and non-evolved stars, which favours 
the primordial scenario. 

In the following, we explore the first two scenarios us- 
ing our fast-rotating models. The abundance pattern ob- 
served at the surface of CEMP stars seems to be a mixture 
of hydrogen and helium burning products. Since rotation 
allows these products to coexist in the outer layers of stars 
(both in massive and intermediate mass stars) , this seems 
a useful direction for our research. 

6.3. Comparison with wind composition 

Let us first see if the CEMP stars could be formed from 
material made up of massive star winds, or at least heavily 
enriched by winds of massive stars. At first sight, such a 
model might appear quite unrealistic, since the period of 
strong stellar winds is rapidly followed by the supernova 
ejection, which would add the ejecta of the supernova itself 
to the winds. However, for massive stars at the end of 
their nuclear lifetime, a black hole, that swallows the whole 
final mass, might be produced. In that case, the massive 
star would contribute to the local chemical enrichment of 
the interstellar medium only through its winds. Let us 
suppose that such a situation has occurred and that the 
small halo star we observe today was formed from the 
shock induced by the stellar winds with the interstellar 
material. What would its chemical composition be? Its 
iron content would be the same as the iron abundance 
of the massive star. Indeed, the iron abundance of the 
interstellar medium would have no time to change much 
in the brief massive star lifetime, and the massive star 
wind ejecta are neither depleted nor enriched in iron. 

The abundances of the other elements in the stellar 
winds for our two rotating 60 Mq at Z = 10~* and 10~^ 
are shown in Fig. (S] The ordinate [X/Fe] is given by the 
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Fig. 8. Chemical composition of the wind of rotating 
60 M0 models (solid circles and squares). The hatched 
areas correspond to the range of values measured at the 
surface of giant CEMP stars: HE 0107-5240, [Fe/H]~ - 
5.3 (Christlieb & al. 2004); CS 22949-037, [Fe/H]~ -4.0 
(Norris & al. 2001, Depagn e fc a l. lM)^ : CS 29498-043, 
[Fe/H]~ -3.5 (Aoki & al. EUUljl . The empty triangles 
(Plez & Cohen [Fe/H]~ -4.0) and stars (Frebel et 
al. I2005L [Fe/H]~ —5.4, only an upper limit is given for 
[0/Fe]) correspond to non-evolved CEMP stars (see text). 

following expression: 

where X is the mass fraction of element X in the wind 
ejecta, and Xq in the Sun. Similarly, the symbols -'^(Fe) 
and X(Fe)o refer to the mass fraction of ^^Fe in the 
wind material or in the Sun. Here we suppose that 
log(X(Fe)/X(Fe)Q) ^ [Fe/H], since the mass fraction of 
hydrogen remains approximately constant, whatever the 
metallicity between Z =10~* and 0.02. The values of 
[Fe/H] are those corresponding to the initial metallicity 
of the models (for Z = 10"^ one has log(X(Fe)/X(H)) = 
—9.38, see Table^). The solar abundances are those chosen 
by Christlieb et al (5001) and Bessell et al. in their 

analysis of the star HE 0107-5240, and they correspond 
to the solar abundances obtained recently by Asplund 
et al. (EOnS- In particular, log{X (Fe) / X (R)) q = -2.80, 
thus [Fe/H]=-6.6 at Z = 10"*, and [Fc/H]=-3.6 at Z = 
10-^ 

From Fig. |S1 we see that the winds are strongly en- 
riched in CNO elements. The model at Z — 10~^, com- 
puted with an extended nuclear reaction network, allows 
us to look at the abundances in the winds of heavier el- 
ements. The wind material is also somewhat enriched in 
Na and Al. Before comparing with the observations, let us 
first note: 

— 1) The more metal-poor model is shifted toward higher 
values compared to the metal rich one. If, in both mod- 
els, the mass fraction of element X were the same in 
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the wind ejecta, then one would expect a shift by 3 dex 
when the iron content of the model goes from [Fe/H]=- 
3.6 to [Fe/H]=-6.6. The actual shift is approximately 
3.6 dex, slightly more than the iron content difference 
of 3 dex between the two models. The additional 0.6 
dex comes from the more efficient mixing in the metal- 
poorer models. 

— 2) In the frame of the hypotheses made here, i.e. a 
halo star made from the wind ejecta that triggered its 
formation, we should compare the wind composition 
from a massive star model with the same initial iron 
content as in the halo star we considered. The range 
of iron contents in the models, [Fe/H] equal from - 
6.6 to -3.6, covers the range of iron content of the 
CEMP plotted in Fig. il whose [Fe/H] are between 
-5.4 and -3.5. However, the [Fe/H] = -6.6 model is well 
below the lower bound of the observed [Fe/H], making 
this model less interesting for comparisons with the 
presently available observations. In that respect the 
[Fe/H] = -3.6 model, which has an iron content that is 
comparable to the iron-richest stars observed, is more 
interesting. 

— 3) Any dilution with some amount of interstellar mate- 
rial would lower the abundance of the element X with- 
out changing the mass fraction of iron. In that case 
the values plotted in Fig. |S1 are shifted to lower val- 
ues, but the relative abundances of the elements will 
not change (as far as the main source of the elements 
considered are the wind ejecta). 

Keeping in mind these three comments, it appears that 
what has to be compared with the observations are more 
the relative abundances between the CNO elements than 
the actual values of the [X/Fe] ratios, which will depend 
on the initial metallicity of the model considered, as well 
as on the dilution factor. 

From Fig. |S1 and Tables 13 & El one sees that, for the 
two metallicities considered here, the wind material of ro- 
tating models is characterised by N/C and N/0 ratios 
between 1 and 40 and by ^^C/^^C ratios around 4-5. 
These values are compatible with the ratios observed at 
the surface of CS 22949-037 (Depagne et al.ESEI): N/C 
- 3 and 12 C / 13 C -4. The observed value for N/0 
(~0.2) is smaller than the range of theoretical values, but 
greater than the solar ratio 0.03). Thus the observed 
N/0 ratio also bears the mark of some CNO processing, 
although slightly less pronounced than in our stellar wind 
models. 

On the whole, a stellar wind origin for the material 
composing this star does not appear out of order in view 
of these comparisons, all the more so if one considers that, 
in the present comparison, there is no fine tuning of some 
parameters in order to obtain the best agreement possible. 
The theoretical results are directly compared to the ob- 
servations. Moreover only a small subset of possible initial 
conditions has been explored. 

Other CEMP stars present, however, lower values for 
the N/C and N/0 ratios and higher values for the ^^C/^^C 



ratio. For these cases, it appears that the winds of our ro- 
tating 60 M0 models appear to be too strongly CNO- 
processed (too high N/C and N/0 ratios and too low 
i2Qyi3(^ ratios). Better agreement would be obtained if 
the observed abundances also stem from material coming 
from the CO-core and ejected either by strong late stellar 
winds or in a supernova explosion. 

6.4. Expression of abundance ratios in total ejecta 
( winds and supernova ) 

To find the initial chemical composition of stars that 
would form from a mixture of wind and supernova ejecta 
with interstellar medium material, let us define X^j as the 
mass fraction of element X in the ejecta (wind and su- 
pernova). This quantity can be obtained from the stellar 
models and computed according to the expression below: 



"T-wind + msN 



where X^ind and Xsn are the mass fractions of element X 
in the wind, in the supernova ejecta. Here mwind and rnsN 
are the masses ejected by the stellar winds and at the time 
of the supernova explosion. To obtain the mass of the rem- 
nants, we adopted the relation obtained by Arnett H1991|l 
between the masses of the remnant and the CO core. This 
method is the same as the one adopted by Maeder l(TU^ . 

The total mass ejected by the star, moj = Tiwind+^SN, 
is mixed with some amount of interstellar material toism • 
The mass fraction of element X in the material composed 
from the ejecta mixed with the interstellar medium will 
be 



-t- 1 



where Xi^i is the mass fraction of element X in the in- 
terstellar medium. In our case the interstellar medium is 
very metal poor so that one can consider Xi^i ^ for the 
heavy elements synthesised in great quantities by the star 
(note that this cannot be done for nitrogen ejected by the 
non-rotating 60 Mq stellar model). We also suppose that 
mcj ^ misM and thus X = (Arcjmcj)/misM- Using these 
expressions, one can write 



[Fe/H] = log 



( XjFe), 



■)+logf^), 

/ V^"-ISM/ 



assuming, as we did above, that X(R)Q/X(R)cj 
Values of [X/Fe] are obtained using the expression 



[X/Fe] 

'Xr. 



[X/H] 



log ( ^ ) - log 



[Fe/H] 
X{FcU 



One needs to have an estimate for both rricj/misM (the di- 
lution factor) and for the mass fraction of iron in the ejecta 
X(Fe)cj. A precise quantitative determination of X(Fe)oj 
and mej/misM from theory is quite difficult. For instance. 
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for a given initial mass, the quantity of iron ejected by 
the supernova can vary by orders of magnitudes depend- 
ing on the mass cut, the energetics of the supernova, and 
the geometry of the expfosion (see e.g. Maeda & Nomoto 
I2003|l . On the other hand, the dilution factor will depend 
on the energetics of the supernova, among other parame- 
ters. In the absence of any very precise guidelines, we de- 
termined the two unknown quantities, the dilution factor, 
and the mass of ejected iron, requiring that the mixture 
will have [Fe/H]= -5.4 and [0/Fe]= +3.5. The first value 
corresponds to the values observed at the surface of the 
star HE 1327-2326 (Frebel et al. ^3^, and the second 
value is below the upper limit of [0/Fe] (< 4.0) found for 
this star. Doing so, one can write. 



[X/Fe] 



log 



[X/0] + [0/Fe] = 

-(it)-- 



where ^(O) is the mass fraction of oxygen. The mass frac- 
tion of ejected iron can be estimated from 



[0/Fe] = log 



^ U(Fe) 



and the dilution factor can be obtained from 



[Fe/H] = log 



fX{Fe 



\X{Fc), 



log 



3.5, 



-5.4. 



6.5. Results from the "wind plus supernova ejecta" 
model 

Using the above formulae, let us now discuss what can be 
expected for the chemical composition of a very metal- 
poor star formed from such a mixture. This includes the 
CNO elements for which the present models can give con- 
sistent estimates and the case of our models at Z = 10~^ 
([Fe/H]=-6.6), which are the models that are compatible 
with the requirement that the mixture of ejecta and ISM 
material must have an [Fe/II]= -5.4. Obviously our second 
series of models at Z ^ 10"^ ([Fe/II]= -3.6) does not fit 
this requirement. Imposing [0/Fe]=3.5 and [Fe/H]= -5.4 
implies masses of iron that are being ejected on the order 
of 1 X 10"'^ Mq and mixed with a mass of interstellar 
medium of about 2 x 10^ M0. The mass of ejected iron 
(actually in the form of ^^Ni) is very small compared to 
the classical values of 0.07-0.10 Mq. On the other hand, 
this quantity can be very small if a large part of the mass 
falls back onto the remnant (Umeda & Nomoto |2503 ) . 
The mass of interstellar gas collected by the shock wave 
can be related to the explosion energy E'exp through (see 
Shigeyama & Tsujimoto[ 



MisM = 5.1 X lO'^M© 



cxp 



lO^^erg 



A mass of 2 x 10^ 
equal to 4 x 10^^ 
energies released by supernova explosion. Thus imposing 



Mq would correspond to an energy 
erg, i.e., a value well in the range of 




element charge number 

Fig. 9. Chemical composition of the ejecta (wind and su- 
pernova) of 60 Mq models: solid circles and triangles cor- 
respond to models at Z = 10^^ ([Fe/H]=-6.6) with and 
without rotation. The [N/Fe] ratio for the non-rotating 
model is equal to 0, i.e., no N-enrichment is expected. The 
hatched areas (moving from top right down to the left) 
correspond to the range of values measured at the sur- 
face of giant CEMP stars (same stars as in Fig. ^ . The 
empty triangles (Plez & Cohen 2Qf^ and stars (Frebel 
et al. 200S1 only an upper limit is given for [0/Fe]) corre- 
spond to non-evolved CEMP stars (see text). The hatched 
areas (L to R from top) show the range of values mea- 
sured at the surface of normal halo giant stars by Cayrel 
et al. H2004|l and Spite et al. Ij2005l unmixed sample only, 
see text). 



[0/Fe]=3.5 and [Fe/H]=-5.4 does not imply unrealistic 
values for the mass of iron that is ejected and for the 
mass of interstellar medium swept up by the shock wave 
of the supernova explosion. 

The theoretical ratios for the CNO elements are shown 
in Fig. ini and compared with the ratios observed at the 
surface of CEMP stars and of normal giant halo stars by 
Cayrel et al. 1)2004(1 and Spite et al. 1)2005(1 . For oxygen, 
the value of 3.5 is obtained by construction, so it does not 
provide any constraint; however, see the previous para- 
graph. More interesting, of course, are the carbon and ni- 
trogen abundances. One sees that both the non-rotating 
and rotating models might account for some level of C- 
enrichment that is compatible with the range of values 
observed at the surface of CEMP stars. However, only 
the rotating models produce N-rich material at this level. 
Figure 01 also shows that the predicted value for the N/C 
and N/0 ratios from the rotating model appear to agree 
with the observed values of these ratios at the surface of 
CEMP stars. Thus as expected, the addition of material 
from the CO core (here ejected at the time of the super- 
nova explosion) to the wind material (mainly enriched in 
CNO-processed material), reduces the N/C and N/0 ra- 
tios. 
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The theoretical values in the wind-plus-supernova 
ejecta model for the ratio ^^C/^^C are between 100 and 
4000 for the rotating models (see Tables El & El- The value 
predicted by the non-rotating model is much higher, on 
the order of 10^. Compared to the observed values, which 
are between 4 and 60, the value of the non-rotating model 
is in excess by at least seven orders of magnitude. The sit- 
uation is much more favourable for the rotating models. 
In this last case, the predicted values are still somewhat 
too high, but by much lower factors. 

Other proportions probably exist between wind and 
supernova mass ejecta than those considered here, which 
would provide a better fit to the observed surface abun- 
dances of CEMP stars. Also more models need to be calcu- 
lated for exploring the set of initial parameters leading to 
a good correspondence between theory and observations. 
From the large range of results obtained by different initial 
conditions, there is little doubt, that such a given set of 
parameters exists. The ^^C/^'^C ratio appears extremely 
sensitive to input parameters, so it may be a powerful tool 
for a closer identification of the exact nucleosynthetic site. 

As can be seen from Fig.O the abundances observed at 
the surface of the normal giant stars by Cayrel et al. (,2004 ) 
and Spite et al. H2()()5|l are not far from solar ratios, and 
are well below the range of values observed at the surface 
of CEMP stars. Only the subset of stars qualified as un- 
mixed by Spite et al. H2005|l . i.e., presenting no evidence 
of C to N conversion, has been plotted here. Probably 
these stars are formed from a reservoir of matter made 
up of the ejecta of different initial mass stars, convolved 
with a proper distribution of the initial rotation velocities, 
while the C-rich stars require some special circumstances 
involving a few or maybe only one nucleosynthetic event. 

6.6. Chemical composition of the envelopes of E-AGB 
stars 

One can wonder whether intermediate mass stars could 
also play a role in explaining the peculiar abundance pat- 
tern of CEMP stars. For instance, Suda et al. (|2.004) pro- 
pose that the small halo star, observed today as a CEMP 
star, was the secondary in a close binary system. The sec- 
ondary might have accreted matter from its evolved com- 
panion, an AGB star, and might have thus acquired at 
least part of its peculiar surface abundance pattern. 

The physical conditions encountered in the advanced 
phases of an intermediate mass star are not so different 
from the one realised in massive stars. Thus the same nu- 
clear reaction chains can occur and lead to similar nucle- 
osynthetic products. Also the lifetimes of massive stars 
(on the order of a few million years) are not very differ- 
ent from the lifetimes of the most massive intermediate 
mass stars; typically a 7 Mq has a lifetime on the order 
of 40 Myr, only an order of magnitude higher than a 60 
M© star. Moreover, the observation of s-process element 
overabundances at the surface of some CEMP stars also 
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Fig. 10. Chemical composition of the envelopes of E-AGB 
stars compared to abundances observed at the surface of 
CEMP stars (hatched areas). The continuous line shows 
the case of a 7 M© at Z = lO"'"^ ([Fe/H]=-3.6) with Vini = 
800 km s'-^. The vertical lines (shown as "error bars") in- 
dicate the ranges of values for CNO elements in the stellar 
models of Meynet & Maeder |2002) (models with initial 
masses between 2 and 7 M© at Z = 10"^). The thick 
and thin lines correspond to rotating (v-mi — 300 km s^^) 
and non-rotating models. The empty triangles (Plez & 
Cohen and stars (Frebel et al. IM)^ only an up- 

per limit is given for [0/Fe]) correspond to non-evolved 
CEMP stars. 



point toward a possible Asymptotic Giant Branch (AGB) 
star origin^ for the material composing the CEMP stars. 

To explore this scenario, we computed a 7 M© with 
Vini = 800 km s^"'^ aX Z — 10^^, and with the same physi- 
cal ingredients as the 60 M© stellar models of the present 
paper. In contrast to the 60 M© models, the 7 M© stellar 
model loses little mass during the core H- and He-burning 
phase, so that the star has still nearly its whole original 
mass at the early asymptotic giant branch stage (the ac- 
tual mass at this stage is 6.988 M©). This is because the 
star never reaches the break-up limit during the MS phase; 
and, due to rotational mixing and dredge-up, the metal- 
licity enhancement at the surface only occurs very late, 
when the star evolves toward the red part of the HR dia- 
gram after the end of the core He-burning phase. At this 
stage, the outer envelope of the star is enriched in primary 
CNO elements, and the surface metallicity reaches about 
1000 times the initial metallicity. If such a star is in a close 
binary system, there is good chance that mass transfer oc- 
curs during this phase of expansion of the outer layers. In 
that case, the secondary may accrete part of the envelope 
of the E-AGB star. 



^ Note that massive stars also produce s-process elements. 
The massive star s-process elements have low atomic mass 
number (A inferior to about 90) and are known under the name 
of the weak component of the s-process. 
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From the 7 Mq stellar model, we can estimate the 
chemical composition of the envelope at the beginning of 
the thermal pulse AGB phase. Here we call envelope all 
the material above the CO-core. The result is shown in 
Fig. [TUI (continuous line with solid circles). We also plot- 
ted the values obtained from the models of Meynet & 
Maeder H2()()2(l for initial masses between 2 and 7 M© at 
Z = 10"^ and with Uini = and 300 km s~^ 

Before discussing the comparisons with observations, 
let us make two remarks: 1) as was the case for the the- 
oretical predictions of the massive star winds, the values 
given here have not been adjusted to fit the observed val- 
ues but result from the basic physics of the models; 2) the 
initial metallicity of our AGB models ([Fe/H]=-3.6) is at 
the higher range of values of the metallicities observed for 
the CEMP stars. However, based on the results from our 
massive star models at [Fe/H]=-6.6 and -3.6 (see Fig.lHJ, 
we see that the overall pattern of the abundances will 
probably remain quite similar for a lower initial metallic- 
ity. Only a shift toward higher values along the vertical 
axis is expected when the initial metallicity of the model 
is decreased. 

Looking at Fig. EH one can note the three following 
points: 

1. The envelope of rotating intermediate mass stars 
presents a chemical composition in carbon, nitrogen, 
and oxygen that agrees well with the observed value 
at the surface of CEMP stars. In particular, compared 
to the wind material of massive stars (see Fig.lHJ, the 
N/C ratios and N/0 ratios are in better agreement. 
The non-rotating models cannot account for the high 
overabundances in nitrogen and oxygen. 

2. The ^^C/^^C ratios in our rotating models are between 
19 and 2500, with the lowest values corresponding to 
the most massive intermediate-mass star models. The 
non-rotating models give values between 3 x 10^ and 2 
X 10^. Again here, rotating models agree much better 
with the observed values, although very low ^'^Gj^^C 
values (on the order of 4-5, as observed e.g. at the 
surface of the dwarf halo star G77 61, see Plez and 
Cohen [2003 seem to be reproduced only by massive 
star models (wind material only). 

3. For sodium and aluminum, the ratios predicted by our 
7 M© model with Uini — 800 km s^^ fit the observed 
values well. In the case of magnesium, good agreement 
is also obtained. 

Thus we see that the envelopes of AGB stellar models 
with rotation show a very similar chemical composition 
to the one observed at the surface of CEMP stars. It is, 
however, still difficult to say that rotating intermediate 
mass star models are better than rotating massive star 
models in this respect. Probably, some CEMP stars are 
formed from massive star ejecta and others from AGB 
star envelopes. Interestingly at this stage, some possible 
ways to distinguish between massive star wind material 
and AGB envelopes do appear. Indeed, we just saw above 
that massive star wind material is characterised by a very 



low ^^C/^^C ratio, while intermediate mass stars seem to 
present higher values for this ratio. The AGB envelopes 
would also present very high overabundances of "'^^O, ^^0, 
^^F, and ^^Ne, while wind of massive rotating stars present 
a weaker overabundance of ^^O and depletion of ^*0, ^^F, 
and ^^Ne. As discussed in Frebel et al. I)2()()5|l . the ratio 
of heavy elements, such as the strontium-to-barium ratio, 
can also give clues to the origin of the material from which 
the star formed. In the case of HE 1327-2326, Frebel et 
al. (2005) give a lower hmit of [Sr/Ba] > -0.4, which 
suggests that strontium was not produced in the main s- 
process occurring in AGB stars, thus leaving the massive 
star hypothesis as the best option, in agreement with the 
result from i^C/iSC in G77-61 (Plez & Cohen and 
CS 22949-037 (Depagne et al. 1200211 . 

7. Conclusion 

We have proposed a new scenario for the evolution of very 
metal-poor massive stars. This scenario requires no new 
physical processes, as it is based on models that have been 
extensively compared to observations of stars at solar com- 
position and in the LMC and SMC. The changes with 
respect to classical scenarios are twofold and are both in- 
duced by fast rotation: first, rotational mixing deeply af- 
fects the chemical composition of the material ejected by 
the massive stars; second, rotation significantly enhances 
the mass lost by stellar winds. The mass loss rates are 
increased mainly because the mixing process is so strong 
that the surface metallicity is enhanced by several orders 
of magnitude. This leads to strong radiative winds dur- 
ing the evolution in the red part of the HR diagram. The 
strongest mass loss occurs at the very end of the core 
He-burning phase. The proposed scenario may allow very 
massive stars to avoid the pair instability. 

We show that material ejected by rotating models has 
chemical compositions that show close similarities to the 
peculiar abundance pattern observed at the surface of 
CEMP stars. We explored the three possibilities of CEMP 
stars made of: 1) massive star wind material, 2) total mas- 
sive star ejecta (wind plus supernova ejecta), 3) and ma- 
terial from E-AGB star envelopes. Interestingly, from the 
models computed here, one can order these three possibil- 
ities according to the degree of richness in CNO processed 
material. From the richest to the poorest, one has the 
wind material, the E-AGB envelope, and the total ejecta 
of massive stars. The imprints on the abundance pattern 
of CEMP stars are thus not the same, depending on which 
source is involved. There is good hope that in the future, 
it will be possible to distinguish them. 

Other interesting questions will be explored in the fu- 
ture with these rotating metal-poor models. Among them 
let us briefly mention: 

— What is the enrichment in new synthesized helium by 
the first stellar generations? This is a fundamental 
question already asked long ago by Hoyle & Tayler 
H1964|l . A precise knowledge of the helium enrichment 
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caused by the first massive stars (Carr et al. 119841 
Marigo et al. I2UU3|I is important in order to correctly 
deduce the value of the cosmological helium from the 
observed abundance of helium in low metallicity re- 
gions (see e.g. Salvaterra & Ferrara l2003|l . Production 
of helium by the first massive stars may also affect the 
initial helium content of stars in globular clusters. If 
the initial helium content of stars in globular clusters 
is increased by 0.02 in mass fraction, the stellar mod- 
els will provide ages for the globular clusters that are 
lower by roughly 15%, i.e., 2 Gyr starting from an age 
of 13 Gyr (Shi 1199.51 see also the interesting discussion 
in Marigo et al. 2003|. In the case of our rotating 60 
M0 at Z = 10~^, 22% of the initial mass is ejected in 
the form of new synthesised helium by stellar winds. 
Thus the models presented here will certainly lead to 
new views on the question of the helium enrichment at 
very low metallicity, provided of course, they are rep- 
resentative of the evolution of the majority of massive 
stars at very low Z. 

— What are the .sources of primary nitrogen in very 
metal-poor halo stars? Primary is produced in 
large quantities in our rotating models. In our previous 
work on this subject (Meynet & Maeder l^002|l . we dis- 
cussed the yields from stellar models a.t Z = 10~^ with 
Viai = 300 km s~^. Such an initial velocity corresponds 
to a ratio Wini/i^crit of only 0.25. This value is lower 
than the value of ^^0.35 reached by solar metallicity 
models with Uini — 300 km s~^. With such a low ini- 
tial ratio of Uini/t^crit, wc found that the main sources 
of primary nitrogen were intermediate mass stars with 
initial masses around about 3 Mq. However, as already 
shown in Meynet & Maeder ((2002|1 . the yield in ^''N 
increases rapidly when the initial velocity increases. 
As a numerical example, the yield in primary nitro- 
gen for the Z = 10~^, 60 Mq with Uini/fcrit equal to 
0.25 was 7 xlO^* Mq, while the corresponding model 
with Wini/wcrit cqual to 0.65 produces a yield of ~2 
xlO~^; i.e., it increased by a factor of nearly 300 ! 
Interestingly, these high yields of primary ^^N from 
short-lived massive stars seem to be required for ex- 
plaining the high N/0 ratio observed in metal-poor 
halo stars (Chiappini et al. I2005|l . Note that massive 
intermediate mass stars, whose lifetimes would be only 
an order of magnitude higher than those of the most 
massive stars, could also be invoked to explain the high 
N/0 ratio observed in very metal-poor stars. The age- 
metallicity relation is not precise enough to allow us 
to distinguish between the two. 

— How does rotation affect the yields of extremely metal- 
poor stars? Other elements such as ^'^C, ^^O, ""^^O, and 
^^Ne are also produced in much greater quantities in 
the rotating models. Computations are now in progress 
for extending the range of initial parameters explored 
and to study the impact of such models on the produc- 
tion of these isotopes, as well as on s-process elements. 



We think that the fact that stars rotate, and may even 
rotate fast especially at low metallicity, has to be taken 
into account for obtaining more realistic models of the 
extremely metal-poor stars that formed in the early life 
of the Universe. 
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